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\ The magnetization process of the S—l/2 antiferromagnetic spin ladder at T=0 is studied 

+^ by the exact diagonalization of finite clusters and size scaling analyses. It is found that a 

C/5 ' 

magnetization plateau appears at half the saturation value (m = 1/2) in the presence of a 
sufficiently large next-nearest-neighbor exchange interaction to yield the frustration, when the 
rung coupling is larger than the leg one. The phase diagram at m = 1/2 is given by the 
analysis based on the conformal invariance. The magnetization curves are also presented in 
q . several cases. 
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£NJ , The spin ladder has attracted a lot of interest in the fields of the magnetism and high temperature 
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X 



superconductivity. Many studies have revealed that the low-temperature properties of a system are 



characterized by the spin gapS'l'S due to large quantum fluctuations. Recently several theoretical 
approachesBIl* suggested the possibility of another spin gap induced by the external magnetic field 
generally in one-dimensional antiferromagnets. The spin gap is then observed as a plateau in the 
magnetization curve. The condition of the appearance of the plateau was rigorously given@) as 

^ ' Q(S — m) = integer, (1) 



where Q is the spatial period of the ground state measured by the unit cell. S and m are the total 
spin and the magnetization per unit cell, respectively. In fact it was theoretically revealed that 
the plateau could appear in some generalized spin ladders; three- leg,& dimerizedB and zigzagi'i) 
ladders. The last one is also equivalent to the bond-alternating chain with the next-nearest neighbor 
exchange interaction which is one of the realistic models of CuGeC>3 . We consider another scenario 
of the plateau in the spin ladder, which is due to the symmetric next-nearest neighbor exchange 
interaction. This implies that both diagonal exchanges are equivalently taken into account at every 
plaquette of the ladder. In contrast to the zigzag ladder, this model has reflection symmetry along 
the rung, realized in the standard spin ladder. A strong coupling approach!) showed that the plateau 



l 



also appears in the spin ladder with the symmetric next-nearest neighbor exchange interaction, in 
the large rung coupling limit. In this paper, the magnetization process of the 5 = 1/2 frustrated 
antiferromagnetizatic spin ladder is studied quantitatively using the exact diagonalization of finite 
clusters and size scaling techniques, to verify the realization of the plateau and give the phase 
diagram at m = 1/2. 

We consider the S=l/2 antiferromagnetic spin ladder with the next-nearest-neighbor exchange 
interaction in a magnetic field described by the Heisenberg Hamiltonian 



H — H + Hz 

L 

Hq = Jl / j(Sl,j • Si^+i + S2,i • S2,i+l) 



(2) 



+ J± E( S l,i • S 2,i) 

i 

L 

+aJi / j(Si,j • S2,i+i + S2,j • Si,i+iy 



H z = -HY J {Si i + S [ 



2,i)i 



(3) 
(4) 



under the periodic boundary condition, where J±, J± and aJ\ denotes the coupling constants of 
the leg, rung and next-nearest neighbor (diagonal) exchange interactions, respectively. (Fig. l)We 




Fig. 1. Spin ladder with next-nearest neighbor exchange interactions along the diagonals 



put Jj_=l in the following. TLz is the Zeeman term where H denotes the magnetic field along the 
z-axis and the energy is specified by the eigenvalue M of the conserved quantity J2i (Sf % + S% {)■ We 
calculate, by the Lanczos algorithm, the lowest eigenvalue E(L, M) of the Hamiltonian Hq in the 
subspace of constant M=J2i (Sfi + 5| J- (L is even and L < 12.) The macroscopic magnetization is 
represented by m = M/L. The magnetic excitation gap giving 5M = ±1 for the total Hamiltonian 
H is given by 

A± = E(L,M±1) - E(L,M)tH. (5) 
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The width of the plateau in the magnetization curve given by 



A = A + + A_ 

= E(L,M + 1) +E(L,M - 1) - 2E(L,M), (6) 

is a useful order parameter .0) When the system has a plateau, A should converge to a finite value 
A(oo)(t^0) in the infinite length limit and the scaled gap LA diverges in proportion to L. In 
contrast, when the system has no plateau, A should decay as 1/L(L— >oo) and the scaled gap LA 
is independent of L. 

First, we investigate the possibility of the appearance of the magnetization plateau at m=M/L = 
1/2. The scaled gap LA of finite systems (L = 8~12) with fixed J\ (=0.4) is plotted against a in 
Fig. 2(a). In Fig. 2(a), it is apparent that LA increases with increasing L for a > 0.5, while it is 
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Fig. 2. (a)Scaled gap LA vs a (b)Conformal anomaly c and exponent r\ vs a(Ji = 0.4) 



almost independent of L around the region a < 0.5. This implies that the plateau would appear 
for large a, while the system is still gapless for small a. The phase boundary is expected to be 
characterized by the Kosterlitz-Thouless (KT) transition© Let us introduce the critical exponent 
T] defined by the spin correlation function (S^S^) ~ (— l) r r~ v . Next we estimate the conformal 
anomaly c in the conformal field theory (CFT)S According to the CFT, if the system is massless 
at m, the size dependence of the energy has the asymptotic form 



(L - oo), (7) 
where e(m) is the energy of the bulk system, c is the conformal anomaly, and v s is the sound 
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velocity. v s is estimated by 



±-[E kl (L,M)-E(L,M)]+o(± 



(8) 



Where E kl (L, M) is the lowest energy with the momentum k± = 2tt/L. CFT predicts the relation 
A ~ 2irv s r]/L (L — > oo), the exponent 77 can be estimated by the form© 



A 



■' E kl (L,M)-E(L,M) +0 {^J- (9) 
It is expected that the conformal anomaly holds c = 1 in the gapless region. In Fig. 2(b), c 
begins to deviate from the value of 1 around a~0.5, indicating the appearance of a plateau there. 
Around the same a value, r\ = 1 happens to occur. An analysis based on the two-dimensional 
sine-Gordon modelQ) suggested that the exponent ij should be 1/n 2 (n is an integer) at the KT 
phase boundary. In ref.14, the isotropic point of the S = 1/2 XX Z chain was also mentioned as 
a good example of the KT transition with rj = 1. We can therefore use the condition 77 = 1 to 
determine the phase boundary below. The boundaries obtained for L = 10 and 12 are shown in Fig. 
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Fig. 3. Phase diagram in the a-Ji plane at m = 1/2 



3. The L dependence is so small that the calculated boundary is believed to be very close to the 
bulk one. Besides the KT line, we also found the first-order phase boundary, which was detected 
as the energy level crossing in the ground state at m = 1/2. Since it is almost independent of L, 
we show the boundary in Fig. 3 only for L = 12. The scaled gap analysis similar to Fig. 2(a) 
revealed that another gapless phase (gapless2) with c = 1 exists above the boundary of the plateau 
for larger a. Thus the phase diagram for m = 1/2 consists of the two gapless phases and the 
plateau phase. The first order line corresponds to the boundary between the two gapped phases 
in the nonmagnetic ground state, which was obtained by the density matrix renormalization group 
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analysis.Ej) According to the present numerical diagonalization approach, the boundary is almost 
independent of the magnetization m even in the massless phase. In ref.15, the two massive phases 
were called the singlet and Haldane phases, which correspond to the gaplessl and gapless2 phases 
in the magnetic state, respectively. Schematic ground states in the classical limit of gaplessl and 
gapless2 are shown in Figs. 4 and 5. The present analysis suggests that the plateau phase at 
m = 1/2 can appear only in the magnetization process from the singlet (not Haldane) phase. The 
two gapless phases can be easily explained in the classical limit, where they correspond to two 
different ordered phases depending on whether the rung or diagonal interactions are dominant. 




> H 



Fig. 4. Ground state in the classic limit for 0.5 J± > a Ji (gaplessl) 



In the quantum system, instead of each order, the characteristic spin correlation function exhibits 




> H 

Fig. 5. Ground state in the classic limit for 0.5 J± < aJi(gapless2) 

the power-law decay and the boundary is slightly deviated from the classical one {aJ\ = 0.5Jj_) 
because of the quantum fluctuation. In Fig. 3, a necessary condition for the appearance of the 
plateau in the present scenario is revealed to be that the rung coupling is larger than the leg one 
(Ji < JjJ. We derive the magnetization curve in the bulk system from the energy E(L,M) of 
the finite systems. In the ground state magnetization process, the magnetic field to change M 
into M + 1 (M - 1 into M) for finite systems is given by H + (L,M) = E(L,M + 1) - E(L,M) 
(H-(L,M) = E(L,M) - E(L,M - 1)). We define the thermodynamic limit of H±(L,M) as 
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H±(m). When the system is gapless at m, the magnetic field to realize m in the bulk system is 
obtained by H(m) = H + (m) = H_(m) and the finite size correction is proportional to 1/L. On 
the other hand, for the gapped state, H + {m) and H_{m) are different and correspond to both 
edges of the plateau, respectively (H-(m) < H+(m)). In Fig. 6, H+(L,M) and H-(L,M) for 
m=l/6, 1/4, 1/3, 1/2, 2/3 and 3/4 are plotted versus 1/L for the parameters (Ji,a) = (0.5,0.6) 
in the plateau phase at m = 1/2. The plot reveals that H + (L,M) and H^(L,M) converge to 
the same value H{m) in the infinite length limit except for m = and 1/2. In contrast it clearly 
indicates ^ which implies that the plateau appears at m = 1/2. So we estimate 

H±{\) by Shanks transformation© which is useful for an exponentially converging sequence. The 
transformation applied for a sequence P n is defined as 



PL 



Pn-1 + Pi 



n+1 



2P n 



(10) 



Using it, we also estimate H + (0) (the critical field such that the first spin gap vanishes). Finally, 
we present the magnetization curve obtained by the above size scaling method detailed in ref.10, 
for (Ji,a) =(0.5, 0.6) and (0.5, 0.8) in the plateau phase in Fig. 7. Both curves exhibit a distinct 
plateau at m = 1/2 and its width increases with increasing a, as well as the plateau at m = 0. 




Fig. 6. H+(L,M) and H-(L,M) versus 1/L 



The plateau at m = 1/2 in the present system corresponds to the case when Q = 2, S = 1 and 
m = 1/2 are held in the condition for quantization of the magnetization (1). This implies that the 
translational symmetry is spontaneously broken and the double periodicity of the original unit cell 
is realized in the ground state. It is expected that the frustration due to the next-nearest neighbor 
exchange interaction stabilizes the spin structure where the singlet and triplet bonds along the 
rung are alternately placed in the plateau phase shown in Fig. 8, as was discussed for the bond- 
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Fig. 7. Magnetization curves 




Fig. 8. Ground state in the plateau at m = 1/2, singlet and triplet bonds alternate along the rung 



alternating chain with the next-nearest neighbor exchange interaction .B* If such a site is realized, 
the two energy levels with the momentum k = and tt should be degenerated. The gap between 
the two levels for m = 1/2 is given by 

A n = E 7T (L,L/2)-E(L,L/2). (11) 

The scaled gap LA^ for J± = 1 and J\ = 0.4 (the same parameter as Fig. 2) is plotted versus a in 
Fig. 9. A,,- is revealed to decay faster than 1/L with increasing L, which indicates the degeneracy 
of the two levels in the thermodynamic limit. It justifies realization of the ground state specified 
by Fig. 8 at m = 1/2. The plateau with Q = 2 was also predicted® for the spin ladder with the 
four spin exchange interaction,® but the mechanism would be different from the present case. 

In summary, the exact diagonalization and the size scaling analyses based on CFT suggest that 
the plateau appears at m = 1/2 in the magnetization curve of the spin ladder, if J± > J\ is 
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Fig. 9. Energy gap between the two levels with k — and -k at m = 1/2 



satisfied and the next-nearest neighbor exchange interaction is sufficiently large (about half the leg 
coupling). We anticipate the observation of the plateau in some real materials. The magnetization 
measurement of a rung-dominant spin ladder such as Cu2(C5Hi2N2)2CLi (ref.18) at a high pressure, 
which might enhance the next-nearest neighbor exchange interaction, would also be interesting. 
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